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cn • Abstract 

>, 
00 , We investigate the scale-dependence of /nl in the self-interacting curvaton model. 

We show that the scale-dependence, encoded in the spectral index Uf^^, can be 
l/^ • observable by future cosmic microwave background observations, such as CMBpol, 

[-^ . in a significant part of the parameter space of the model. We point out that together 

^P I with information about the trispectrum 5nL; the self-interacting curvaton model 

parameters could be completely fixed by observations. We also discuss the scale- 
dependence of (7nl and its implications for the curvaton model, arguing that it could 
provide a complementary probe in cases where the theoretical value of nf^^ is below 
p\ ' observational sensitivity. 



1 Introduction 

It is now widely recognized that non-Gaussianity is a useful tool to distinguish between 
various models for the primordial perturbation. Although the power spectrum of the CMB 
temperature fluctuations are well constrained by current observations such as WMAP [1], 
at linear order many mechanisms can give rise to almost the same spectrum of fluctuations. 
Thus, to differentiate between models, one needs to go beyond the linear order. For 
example, the standard slow-roll inflation models predict an almost Gaussian primordial 
perturbation, while other mechanism such as the curvaton [2}B], yield almost identical 
power spectrum but can have a significant non-Gaussian component. 

The deviations from the purely Gaussian perturbation are usually characterized by 
some non-linearity parameters, such as /nl, related to the 3-point function, and tnl and 
(7nl, which are related to the 4-point function. The present observational constraint on /nl 
is —10 < /nl < 74 (95%C.L.) [T|^H If taken at face value, the current limit might actually 
suggest that /nl deviates from zero (although it is still consistent with zero); this might 
spell trouble for the standard single-field slow-roll inflation models, which typically give 
/nl which is proportional to the small slow-roll parameters, or at most /nl ~ C'(l) [TO] . 
If this suggestion turns out to be confirmed by Planck, the curvaton models would be 
natural candidates for the primordial perturbation. 

Indeed, non-Gaussianity in the curvaton model has been the subject of extensive study 
pH - [29] . However, in the majority of these studies the curvaton potential is assumed to 
be of the quadratic form. Since in order to be able to decay, the curvaton has to couple 
to some other fields, it then follows that through loop corrections there will also appear 
some effective curvaton self-couplings, even if they were absent at tree level. Moreover, it 
has been pointed out that even small deviations from the quadratic potential may have 
significant effects for the primordial perturbation if due to the curvaton [T2 | [T5 |l2T| - [28] . In 
particular, it has been shown that for non-quadratic curvaton potentials the predictions 
for /nl and (^nl can be quite different as compared with the purely quadratic case [151 
|2T] . In the presence of non-linearities, it is also possible that the non-Gaussianity of the 
perturbation can be revealed by (^nl (or the 4-point function) rather than by /nl, which by 
virtue of the non-linearites can remain very small for certain choices of the parameters [2Tj . 

Another point worth stressing is that the well-known relation /nl ~ l/'^dcc for the cur- 
vaton model, where rdec is related to the curvaton fraction of the total energy density at 
the time of curvaton decay, is not generic and strictly speaking only true for the quadratic 
potential [2S]- Thus a general, self- interacting curvaton model is a source for new obser- 
vational signatures. However, self-interaction means also more free parameters. Hence it 
would be desirable to have more observational quantities beyond /nl and (7nl that would 
help to constrain the models. 

Indeed, recently it has been pointed out that the scale-dependence of /nl could be 
measured by future observations [SO]- Theoretical aspects of the scale-dependence have 
been investigated in [HTHM] . We note that the scale dependence of equilateral type non- 



"^^ For the constraints on tnl and ^nl, see [5H7] and for forecasts on future constraints see [HUH]- 



Gaussianity is also of interest [361439] . Here we study the scale-dependence of the non- 
Gaussianity parameters in a self-interacting curvaton model. We show that, in some cases, 
the scale-dependence of /nl can be large enough to be detected in future observations. 
We also demonstrate that by using the combined information about /nl and (?nl, we may 
essentially fix the curvaton model parameters. 

The organization of this paper is as follows. In the next section, we summarize the 
formalism for the study of the scale-dependence of /nl, which manifests itself in the 
spectral index Uf^^, in a general curvaton model. Then in Section [3l we closely look 
at n/j^L i^ various cases and discuss its detectability in future observations. We also 
discuss the implications for (yfNL, in particular focusing on the complementary nature of 
the combined information. The final section is devoted to the summary and a discussion 
about the implications of the results. 

2 Formalism and Definitions 

For the purpose of discussing the scale-dependence of non-Gaussianity parameters in the 
curvaton model, let us assume that the potential for the curvaton a can be written as 

where rrirj is the mass of the curvaton, and A and p represent the strength and the power 
of the self-coupling, respectively. To characterize the relative size of the self-interaction 
term in the potential, we define the parameter s as follows: 
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where the subscript * represents that the quantity is evaluated at the time the scale leaves 
the horizon. Thus s is the ratio of two terms in Eq. ([T]) at the horizon exit. Once the 
curvaton begins to oscillate, the energy density of the curvaton can soon be well described 
by the purely quadratic case (i.e., the "nearly" quadratic case discussed in [T5l|2T]), the 
curvature perturbation up to the third order, adopting sudden decay approximation, can 
be written a; 
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"^"^ In the following, we do not eonsider isocurvature fluctuations. For studies of isocurvature fluctuations 
in the curvaton model, see e.g. 



where raec is defined as 
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We assume that, for a given infiationary Hubble parameter H^,, the curvaton decay rate 
can be adjusted such that the correct perturbation amphtude ( ~ 10~^ can be obtained. 



The non-hnearity parameter /nl is then given bjl^ 
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Let us now discuss the scale dependence of /nl in the self-interacting curvaton model. 
The scale dependence of /nl is defined as (see [33] for details) 
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where we regard n/j^^ ^s the spectral index for /nl- From this definition, treating Uf^^ as 
a constant is equivalent to assuming the power-law form: 



/nlOc/c'^^nl. 



(7) 



For the self-interacting curvaton model, the spectral index can be calculated as [33] 
(for a very recent, alternative method of calculation see |34] ) 
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where a prime indicates the derivative with respect to cr,,, i.e., V = d^V/da^. 

In addition to the spectral index Uf^^^ , the running of nj^^^ can also be defined analo- 
gously to the running of the spectral index of the usual spectrum as 
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For the curvaton model with self-interaction, we obtain the following expression for the 
running aj,^^ 
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"^^ For the definition and details of the calculations of the non-linearity parameters in the model, we 
refer the reader to e.g., [2T] . 

"^^ We thank Qing-Guo Huang for pointing out an error in this formula in previous versions of this 
paper. The result given here agrees with the arXiv version of |35l . 



The derivative of the potential can be evaluated once s,p,X, and the mass rricr are 
given. However, in the following discussion, instead of the mass we find it convenient to 
use the parameter 770-0-, defined as the ratio between V" and the Hubble parameter during 
inflation: 

V"(cr*) "^a A , pip- 1) 



Notice that the spectral index for the power spectrum in this model is given by 

n, -l = -2e + 2r/,,, (12) 

where e is the slow-roll parameter for the infiaton sector. Thus the value of rj^^ should be 
small enough to be consistent with the observational constraint: 0.9358 < Ug — l < 0.9921 
(95% C.L.) p. 

With the definition (ITT]) for 770-0- , we can express the higher derivatives of the potential 

as 

V'" _ p{p-l){p-2)s 1 

3H^ ''^^ 2 + p{p-l)s a.' ^ ^ 
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It should be noticed here that, in the self-interacting curvaton model, the value of 
/nl can change its sign when parameters are varied. In fact, when /nl changes its sign, 
the scale-dependence of /nl cannot be well described by the power-law form of Eq. (jTj). 
One could also see this by noticing that the expression of Uf^^^ given by Eq. ([8]) diverges 
when /nl = 0. However, as we will discuss later, when /nl is small enough, it would 
be impossible to measure the scale-dependence of /nl; for small enough /nl , it is not 
an observable. Furthermore, we have checked that the magnitude of the running a/j^^ is 
much smaller than Uf^^^ except in the region where /nl ^ C^(1)- In the following we do 
not consider the unobservable region /nl ^ C'(l); as a consequence, for us Uf^^ is indeed 
the proper way to describe scale-dependence of /nl in the whole range of interest. 

Although in this paper we mainly focus on the scale-dependence of /nl , we can also 
define the scale dependence of ^^nl similarly to the case of /nl as 
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Since in the curvaton model with a self-interaction term, the non-linearity parameter (^nl 
is given by 
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as in the case for nf^^, a similar calculation (consistent with the results of [M]) yields n 
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3 n/j^^ in the curvaton model 

We have scanned the parameter space and evaluated the spectral index n/j^^ numerically. 
In the left panel of Fig. [H we plot the value of Uf^^^ as a function of the self-interaction 
power p . In the figure, we have chosen the parameters as s = 0.05 (red solid line), 0.07 
(green dashed line) and for the rest of this paper we choose 

r/,^ = 0.005. (18) 

Note that the scale dependence of both Uf^^ and n^j^^ is proportional to rj^-a, so they would 
both be larger if we chose a larger rj^^ . The corresponding values of /nl are also shown in 
the right panel. We have chosen the value of rdec such that we obtain /nl = 50 at p = 6 
for each relative self-interaction strength s (therefore the values of rdec are different for 
each case). As discussed in the previous section, when /nl very small, with /nl ^ C){1), 
the power- law form cannot well describe the scale dependence of /nl- To emphasize this 
fact we cut out the regions with /nl < 10. 

According to [SU], future observations such as CMBpol may be able to detect the scale 
dependence of /nl- The observational error has been estimated as (for the fiducial value 

of ^/nl = 0) 

50 1 
Anj^^ ^ 0.05- ^=, (19) 

Jnl yjsky 

where /sky is the fraction of the sky measured. Hence in the parameter space where 
the theoretical estimate nj^^ > An/j^^^ we may be able to constrain the self-interacting 
curvaton model by the scale dependence of /nl- For example, as can be read off from 
figure [21 when p = 6 and s = 0.07, the spectral index is Uf^^^ ^ 0.1, which should be 
observable with future observations. 

To chart the parameter space where rif^j^ is a useful tool for model building, we plot 
in Fig. [2] the contours of n/j^^ ^^^ /nl in the s-rdec plane for the cases of p = 6 and 8. 
From the expression of Eq. (JSj), one can see that while generally /nl is proportional to 
1/rdcc, in the self- interacting curvaton model /nl also depends on s through the non-linear 
evolution of the curvaton field. This is encoded in o"osc and its derivatives. On the other 
hand, when rdec is small, nj^^^^ does not depend on rdec but strongly depends on s. Since 
the scale-dependence of /nl mainly originates from the non-linear evolution of a due to 
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Figure 1: Plots of njj^^^ (left panel) and /nl (right panel) as a function of the power p. 
Here we take s = 0.05 (red line) and 0.07 (green line) with ri„„ = 0.005. The shaded areas 
correspond to the regions where the power-law form /nl oc k^f^i- does not hold. 
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Figure 2: Contours of /nl and Uf^^^ for p = 6 (left panel) and 8 (right panel). The region 
with /nl < 10 is shaded. 
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Figure 3: Regions testable in future observations. 



Shown is the parameter space where 
T9l) . (i.e., region where 
0.05 X (50//nl) is satisfied), for 10 < /nl < 100. 



the theoretical prediction for nj^^ exceeds An/j^^ given in Eq, 



the higher order term in the potential, the value of s controls n/j^^- Notice that, when 
^dec *^ C>(1)) the sign of /nl is determined by the factor 1 + croscO"osc/^osc which appears in 
the first term in Eq. ([5]). This factor depends on s and p, but not on rdec- This is the reason 
why the contours of /nl in Fig. [2] appear to converge at s ~ 0.077 and 0.032 for p = 6 and 
8, respectively, where /nl abruptly switches from a positive value to a negative one. The 
region around this line also corresponds to the breakdown of the power-law description for 
the scale-dependence. By cutting out the region with /nl < 10, which anyway is irrelevant 
from an observational viewpoint, we therefore guarantee that the power-law description is 
valid. 

Since the sensitivity of Uf^^^ in future observations depends on the value of /nl, as 
indicated in Eq. (TT9|) . the self-interacting curvaton model may be tested in certain regions 
of the parameter space. This is depicted in Fig. |3l where the regions with Uf^^ > Auf^^^ are 
shown for the cases p = 4, 6 and 8. Here we conservatively adopt the range 10 < /nl < 100. 
If future observations find a non-zero value for Uf^^^^, the self-interacting curvaton model 
must lie inside the colored regions. On the other hand, if no scale-dependence of /nl will 
be detected, the colored region will be excluded. Planck is expected to have about half 
of the sensitivity to Uf^^^ compared to flT9|) . but large scale structure probes may become 
even more sensitive to nj^^^ than CMBpol. 

4 Constraints implied by the trispectrum: ^nl 
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Figure 4: Plots of nj^^ (right) and (^nl (left) for the cases with {p, s) = 
(8, 0.0214). The value of s is chosen such that the cases of p = 6 and 8 
same values of Uf^^ as a function of /nl- 
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Figure 5: Plot of Ug^^ as a function of the power p for s = 0.07. 



In the previous section, we have shown that n/j^^ can be a useful tool to probe/test 
the self- interacting curvaton model. However, even if we can observe Uf^^^ and /nl, the 
situation remains degenerate: different combinations of the power p, the strength s and 
r^ec can give the same values of /nl and nj^^. To demonstrate this, in the right panel 
of Fig. m we plot Uf^^ as a function of /nl, which is almost equivalent to the plot as a 
function of r^cc, for the cases of p = 6 and 8. To give an example, by fixing s = 0.0622 
and 0.0214 for p = Q and 8, respectively, the lines for p = 6 and 8 almost coincide, as is 
shown in the right panel of Fig. |H However, the predictions for (^nl are in general different. 
This is demonstrated in the left panel of Fig. HI where we show (^nl as a function of /nl 
for the same parameters as in the right panel. This figure indicates that, in principle, it 
is possible break the degeneracy between n/^^ ^^^ /nl by including (^nl in the analysis. 
Thus combining together information about the bispectrum amplitude, the spectral index 
of the bispectrum, and the amplitude of the trispectrum, one may actually fix completely 
the parameters of the simplest self-interacting curvaton model ([T]). We note that tnl in 
this model does not give new information, since it satisfies tnl = (6/nl/5)^ (for more 
discussion on the trispectrum see |l5lll6]), but testing this consistency relation is still 
valuable since it can determine whether it was justified to neglect the the infiaton field 
fluctuations [T9] . 

Furthermore, we can also address the scale-dependence of (j'nl; the appropriate expres- 
sion was already given in Sec. Hi In Fig. Ug^^ is plotted as a function of the power 
p for s = 0.07. It is worth noticing that Ug^^^ becomes larger the smaller p is. This is 
in contrast to the behavior of Uf,^^, which becomes smaller as the value of the power p 
becomes smaller, as is shown in Fig. [H Thus, in future observations it may be possible 
to see the scale-dependence of (^nl even if we cannot detect the scale-dependence of /nl- 



Hence, in addition to (^nl itself, the scale-dependence of ^^nl might also provide a tool to 
probe/test the self-interacting curvaton model. 

5 Conclusion and discussion 

In this paper, we have discussed the scale dependence of the non-linear parameters such 
as /nl and (^nl m a curvaton model with a simple self-interaction term given by the po- 
tential ([T]). The free parameters of the model are the curvaton mass m, the power of the 
self-interaction term p, and its relative strength s. The complete parametrization of the 
curvature perturbation would require also the knowledge of H^, determined by the inflaton 
sector. In principle H^ can be observationally fixed by a detection of primordial gravita- 
tional waves, but their amplitude is likely to be very small in the curvaton scenario. In a 
complete model, the curvaton decay rate might not be a totally independent parameter, 
but here we assume that, for a given H^, the decay rate can be tuned to yield the correct 
perturbation amphtude ( ~ 10~^. In such a situation, as has been argued before, even a 
modest amount of non-linearity in the curvaton potential [T5l[2T] may lead to observable 
consequences. Therefore, although /nl is scale-invariant when the curvaton potential is 
quadratic, in many cases the spectral index Uf^^^ of the bispectrum can be detectable when 
one includes the self-interaction term; this fact is depicted in Fig. [H 

Since the observational sensitivity to the scale-dependence of /nl depends on the mag- 
nitude of /nl, there are regions where the self- interacting curvaton model parameters can 
be pinpointed by the observations of /nl and its scale dependence n/j^^- The testable 
domains in the parameter space are shown in Fig. [3l If no scale-dependence of /nl will be 
found, these regions would be excluded so that for the self-interacting curvaton model, the 
scale-dependence can indeed provide us with clear-cut evidence for or against the model. 

In addition to the scale-dependence of /nl, we also studied the non-linearity parameter 
of the trispectrum, (yfNL, and its scale-dependence n^j^^- -^^ i^ shown in Fig. HI even though 
in principle there is a degeneracy in the space of model parameters in that a given pair of 
values for Uf^^ and /nl is produced by several combinations of (p, s), the degeneracy can 
be broken by s'nl- The spectral index of ^^nl could be used to probe the scale-dependence 
of non-Gaussianity in the self-interacting curvaton model at least for some of the cases 
where nf^^ is not observable. 

So far the studies of non-Gaussianity in the curvaton and other models have focused 
on the non-linearity parameters /nl, tnl and (^nl- Here we have pointed out that the scale 
dependence could also provide a realistic and useful tool to test the models. In particular, 
we have argued that the self-interacting curvaton model parameters could actually be 
completely fixed by the combined observations of /nl, S'nl, and the spectral index rifj^^. 
The fact that the projected observational sensitivity of the latter appears in many cases 
to be much higher than would be required for verification of the theoretical estimates is of 
course encouraging. It also lends some hope for the expectation that the physical origin 
of the primordial perturbation can be determined in the not too distant future. 
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